Precise calculation of the Casimir force between gold surfaces 
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We analyse the result of precise measurement of the Casimir force between bodies covered 
with gold. The values of the parameters used to extrapolate the gold dielectric function to low 
frequencies are very important and discussed in detail. The finite temperature effect is shown to 
exceed considerably the experimental errors. The upper limit on the force is found which is smaller 
than the measured force. Many experimental and theoretical uncertainties were excluded with gold 
. . . covering and we conclude that, possibly, a new force has been detected at smaU separations between 

• bodies. 

o ■ 

O ' 12.20.Ds, 12.20.Fv, 03.70.+k, 14.80.-j 

(N . _ 

^jpj' The Casimir force attracts considerable attention nowadays. At early stage it has been measured between two 

^ I conducting plates with large errors and much more accurately for plate and sphere covered with chromium [^1 . A 

. new wave of interest raised when the force between metallized bodies has been measured with high precision . 

■ Theoretical prediction for the force is not so straightforward as expected. In the original Casimir work [Q the force 

' i^P'(a) = j4o^ has been found as change in the zero point energy between two parallel perfectly conducting plates 
I separated by the distance a. For sphere and plate the expression above has to be modified with the proximity force 
^ '. theorem (PFT) to the following 

■ 360 a3' 

O ; 

OO ', where R is the sphere radius. Real metals are not ideal conductors and there is an important correction to taking 

■ into account their real properties. In the Lifshitz theory which generalizes the Casimir approach, the force is 
I represented in terms of the metal dielectric function e (iui) at imaginary frequencies. This theory was applied for 

description of the experiments in |^,^. Different but equivalent method has been developed in |jl^, where detailed 
calculations were made for Al, Au, and Cu metallization of the bodies using the handbook optical data for these 
Oh, metals. The result for Au is in good agreement with the torsion pendulum experiment [Q. 

' For the atomic force microscope (AFM) experiment Al metallization has been used. This metal is oxidized 
^ I easily and charges can be trapped in the oxide. To prevent influence of these charges on the force, Al was covered 
^ ■ with a thin layer (20 or 8 nm thick) of Au/Pd. The conclusion has been made [3[ that the result of the AFM 
experiment was in agreement with the theory. However, in all calculations [pUGlJl 1 1 the role of Au/Pd layer was 
• • r ignored. Importance of this layer was stressed in ||l^ (sec also ) , where the upper limit on the force has been found 
^ using parameters of single-crystalline materials. This limit was smaller than the measured force in both experiments 
^ ||,| and the discr epancy far exceeded the experimental errors. To exclude uncertainties connected with aluminum 
oxidation, Au/Pd layer and clarify the origin of the discrepancy, it was proposed [T^ j to use Au metallization instead 
5^ , oi Al. Very recently the AFM experiment with Au metallization has been carried out [Q. In this experiment the 
residual potential was reduced to a negligible level so as the surface roughness. The force was measured for smaller 
separations up to 62 nm with uncertainty 3.5 pN . In the previous experiments the photodiods response on the 
cantilever deflection was calibrated assuming that the plate-sphere system on the contact behaves as a rigid body. 
As the new experiment showed this natural assumption was not quite correct and different calibration method was 
proposed [ |l^ . To all appearance the excessive force revealed in the experimental data [ p^ originates mainly from the 
calibration effect. 

The new experiment gives us a solid ground for verification of the theory. Though good agreement with the theory 
is reported ||lj|, this conclusion cannot be considered as final because of the following reasons. Deviations between 
theory and experiment were significant at small separations but the authors averaged them over the whole range of 
separations. An error was admitted in extrapolation of the optical data for Au to low frequencies. At last, the finite 
temperature correction was neglected but it is important at small separations [^,^ . These are the main points we 
will discuss in this paper. 

We start from the general expression for the Casimir force given by Lifshitz |^ modified for the case of sphere-disk 
geometry: 
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Fia) = -^Y1 f dppln[(l-Gi)(l-G2)], (2) 
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Gi = ( £-^e-f?""/M , G2 = ( ^^e-PC"«/^ 



s = «n) - C« = 2nnkT/h. (3) 

e (iCn) is the dielectric function of Am at imaginary frequencies. The prime over the sum sign indicates that the first 
n = term has to be taken with the factor 1/2. Special care needs to treat this term. The formal reason is that 
becomes zero but the integral over p diverges. The physical reason is that this term corresponds to the static limit 
when for metallic bodies we should take e — > 00. The resulting contribution in the force corresponds to the classical 
limit Fci (a) for metals p^ ] 
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F.i{a)^^C{3), (4) 

where C (N) is the zeta-function. If the metal is ideal, this term is canceled by the rest ones in the sum leaving 
only negligible temperature correction. The cancelation is incomplete for real metals and the correction becomes 
important. To evaluate the Casimir force, the integral instead of the sum in (|^) was used in The difference 

between the sum and integral is just the temperature correction AtF which is nearly linear in T. The relative 
value of AtF becomes smaller for smaller separations but its absolute value increases. To have an estimate for 
the correction, we used the Drude approximation for e with the parameters of perfect single-crystalline gold 
ujp = 1.37-10"'^^ s~^, lOt = 3.7-10^^ and found for AtF the values 4 and 14piV at a = 100 and 60 nm, respectively. 
Obviously this correction cannot be ignored in precise calculations. 

The function e (iCn) cannot be measured directly for a given material but can be expressed via imaginary part of 
the dielectric function on the real axis e" (w) with the dispersion relation 
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Actual uncertainty in the force calculations arise from uncertainties in the function e" {uj). Part of the available 
data on e" (w) is shown in Fig . 0. The solid line represents the data taken from ijl^ which were extrapolated to low 
frequencies as discussed in pO| | (see below). The points marked with the open squares and triangles are taken from 
1^ . The circles represent the data from jl^ . The points marked with the crosses are the only data in the figure 
for the evaporated gold film while all the rest are for the bulk material. The data can be separated conventionally in 
two frequency regions with the boundary at the curve minimum loi — 3.2 • 10^^ . In the region oj < uji the electron 
scattering on defects dominates in the material absorptivity. That is why in this range one can find different values 
for e" (w) corresponding to the samples prepared in different conditions. At <C the dielectric function is well 
described by the Drude model 

^2 

£{uj) = e' + ie" = 1-—-^—. (6) 

OJ [Ijj + ILUt) 

The parameters one can find by fitting the data for the real and imaginary parts of e (lu) with (^. It has been done 
for the points marked by squares and triangles (rows 2 and 3 in Table |). Different procedure has been used 
for the solid line |T^. The data in are available for frequencies huj > 0.1 eV. They have been extrapolated to 
lower frequencies according to (|^). The plasma frequency was estimated using its relation with the concentration n of 
free electrons and supposing that every atom produce one electron: tOp = ^ e^n/m*eo- Here e is the electron charge, 
£0 is the free space permittivity, and m* is the effective electron mass which is close but larger than me. Then at 
fixed Wp the data |l^ have been fitted with (^) to find lo^ (row 1 in Table |). The found in this way has actually 
the largest of possible values and, of course, the real material has smaller Wp due to presence of the defects. 

In contrast to low frequencies all the points at w > wi describe practically the same curve. That is because at high 
frequencies the interband absorption dominates and the defects become unimportant. Therefore, one can consider this 
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part of the curve as the universal one (sample independent) . We estimated uncertainty in e [iQ due to a reasonable 
variation of this part of the curve. It is less than 1% for important frequencies C ^ c/2a. The resulting uncertainty 
in the force will be much less than the experimental errors. One can separate contributions in e {iC,) from different 
frequency regions and write e (iC) = 1 + Ei + £2 + £3- The contribution of low frequencies < a; < where the 
optical data are unavailable, can be expressed analytically using (ph in the dispersion relation: 



£1 «) = ^ 



tan — — tan 



£2 and £3 can be found numerically using the data in the ranges loq < lu < uji and w > wi , respectively. To get an idea 
of the contributions of different intervals, we have found for the solid line in Fig. 0£i = 26.6, £2 = 8.3, and £3 = 5.3 
at C = c/2a with the smallest separation a = 63 nm (the first point presented in Let us stress once more that 

£3 does not depend on the defect concentration in the material while Si is quite sensitive to this parameter via actual 
values of ujp and Ur and its contribution dominates in £ («C)- 

Proceeding exactly as in we successfully reproduced £ {iC,) given there in Fig. 2. However, we failed to reproduce 
the same function in ||20|| though the authors claimed to use the same data and the same values for and w^. The 
difference between the curves presented in and is 46% at C = 10^^ and decreases with frequency increase. 
To understand the origin of this difference, consider the limit C ^ when the dielectric function behaves as £ {iQ 
(£oPC)~^- Here p is the static resistivity defined via the Drude parameters as p = lOt/eqUJ^. The parameters ujp 
and lOt used in pO| correspond io p = 3.2 /iil • cm. For unknown reason the function e (iC) presented in [ pO[ at low 
frequencies agrees much better with p = 2.3 n^l ■ cm that has no relation with the data in p^ . 

The Casimir force was evaluated in pc| ] and pO| without the temperature correction using the integral instead of 
the sum in (|^). To check the procedure, we repeated the calculations there at a = 100 nm and found for the reduction 
factor 1] = F/F^ the values 0.547 and 0.559 to be compared with 0.55 jl^ and 0.56 |20|, respectively. The small 
difference in 77 results from different resistivities used to extrapolate the optical data. At a = 63 nm this difference 
is not so harmless: we found for the force F = 467 pN and 477 pN, respectively. The latter value is in agreement 



with the theoretical curve in since the calculations have been done as in 1 20 1 . The force has to be compared with 
the experimental value of 491 pN . The difference is 4 times larger than the experimental errors although the rms 
deviation of the experiment from theory for the whole range of separations is only 3.8 pN p^ . 

Now let us take into consideration the temperature correction which is positive and will improve agreement with 
the experiment. We make calculations using Eq.(|^) with e" (uj) given by the solid line, squares, and triangles in Fig. 
|l|. For the squares there is a gap in the interval 6.3 • 10^"* s^^ < uj < oji which is filled by the linear interpolation 
on a log-log scale. Different ways to fill the gap lead to uncertainty in the force ±1 pN. The results of calculations 
are presented in Table ||. Note that the force given in iQ accidentally coincides with the right one (row 1) because 
the error in p was compensated by the temperature correction. The parameters in row 2 are quite close to those of 
perfect single-crystal. Therefore, the force has to be maximal in this case |lj] due to maximal reflectivity. However, 
we see that it is a bit larger for row 1. The reason is that ujp in ||l^ was not extracted from the optical data but set 
to be maximal by hand. Real material will have smaller cOp that will lead to smaller force. The residual force defined 
as AF (tti) = Fexp (fti) — Ftheorio^i), whcrc Qi are the experimental separations, is shown in Fig. ^ for < 103 nm. 
It has been calculated for the squares in Fig. ^ and can be considered as the minimal residual force. As we can see 
AF exceeds 4 standard deviations in a few points of the closest approach. Parameters of evaporated gold films will 
deviate significantly from those of the perfect single-crystal because the films contain much more defects than the bulk 



material |13|. Therefore, the residual force for real body covering will be definitely larger and will be significant at 
larger separations. One can easily improve the upper limit on the predicted Casimir force by measuring the resistivity 
of the gold films used in the experiment. 

The experiment has been substantially improved in comparison with and many uncertainties have been 
excluded. If there is no an additional systematic error at small separation, the deviation between theory and exper- 
iment has to be considered seriously. In this connection it is very important to reproduce the result in independent 
measurements. Good potential has the micromechanical resonator experiment which already has the necessary 
sensitivity but an additional work has to be done to control the separations ~ 100 nm. 

One can speculate that the residual force is explained by a new Yukawa force mediated by a light scalar boson. In 
this case interaction of two atoms is described by the potential 

he 

Vy (r) = -aiViiVs— exp(-r/A), (7) 
r 

where a is the dimensionless interaction constant, A is the Compton wavelength of a particle responsible for the 
interaction, and A^i^2 is the number of nucleons in atoms of the interacting bodies. With 95% cl the residual force is 
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larger than 17 — 2 • 3.5 = 10 pN . Integrating for the bodies in the experimental configuration one finds the lower 
limit on the interaction constant 

6.27- 10~25e'^/^ flQQnmV' 
" ^ 1 - 1.74e-''A + 0.75e-2VA A ) ' ^ ' 

where a = 63 nm is the minimal separation and h — nm is the gold film thickness. The interaction constant is 
restricted from the He burning stars as a < 1.5 • lO^^s Then the allowed region for the wavelength is A > 33 nni 
or for the scalar boson mass m < 38 eV . 

In |2^ it is stated that the sphere-plate geometry cannot be used for verification of the Casimir effect since a 
rigorous theoretical calculation has never been carried out for this configuration. In this respect we would like to 
note that such a logic does not allow to verify the effect for finite-size plates, either, because the plates always can be 
considered as parts of spherical surfaces. From physical point of view it is obvious that the force for infinite plates 
will be a good approximation il L ^ a, where L is the plate size. The same is true for a sphere and plate if i? 3> a. 
One can discuss the applicability of PFT but there is no any reason to worry (see discussion in and jl^). 

When this work was finished we become aware about the paper [2^ , where the authors claimed that the Schwinger 
prescription [|l^ for the n ~ term in (|^) cannot be applied for nonideal metals. As the result they found the force 
in the classical limit which is 2 times smaller than (^). We think that the statement needs further investigation (see, 
for example, [^), but in any case our conclusion will not change since the Casimir force will be only smaller for the 
smaller n = term. 

In conclusion, we found the upper limit on the Casimir force between gold covered plate and sphere using the 
handbook data for e {lo). The effect of nonzero temperature is shown to be considerably larger of the experimental 
errors. The excessive force found at small separations exceeds 4 standard deviations. The lower limit on the interaction 
constant is found assuming that the residual force is the manifestation of a new Yukawa interaction. 
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FIG. 1. Imaginary part of the dielectric function e"{uj) for gold. The data were collected from different sources (see expla- 
nations in the text). 



FIG. 2. The residual force as a function of separation calculated for the squares in Fig. 
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F{amin) [pN] 
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1.38 


5.38 


3.19 


477 
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1.37 


4.06 


2.44 


474 


3 


1.28 


3.29 


2.27 


459 



TABLE I. The Drude parameters, resistivity, and the Casiniir force at amin = 63 nm for e"(w) represented by the solid line 
(1), squares (2), and triangles (3) in Fig. |l| 
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